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INTRODUCTION 


Litton  Systems,  Inc*,  delivered  inertial  positioning  system  number  one 
(iPS-l)  to  the  U.S.  Defense  Mapping  Agency  in  1975*  The  system  became  oper¬ 
ational  a  year  later.  System  utilization  has  been  extensively  reported  by 
(Harris,  et  al  1976),  (Harris,  1977),  ( Doxey ,  1977),  (Sharp,  1979),  and 
( Blaokmev ,  1981). 

IPS-1  is  the  Litton  Autosurveyor™  with  a  vertical  channel  A-1000 
accelerometer,  as  used  by  the  U.S.  Bureau  of  Land  Management  (BLM),  the 
Canadian  Department  of  Energy,  Mines  and  Resources  (EMR) ,  the  U.S.  Army 
Engineer  Topographic  Laboratories  (ETL),  and  SPAN  International.  The  survey 
data  output  are  raw  and  smoothed  values  for  station  latitude,  longitude  and 
height,  plus  the  relative  free-air  gravity  anomalies  and  output  from  which  the 
relative  deflection  of  the  vertical  components  may  be  computed. 

The  raw  and  smoothed  survey  values  are  stored  in  a  permanent  data  file  of 
the  VAX  11/780  mainframe  computer.  Using  that  file,  certain  post-mission 
processing  is  possible.  For  a  suitably  designed  and  executed  survey  project 
or  area,  it  is  possible  to  improve  the  overall  results" by  a  least-squares 
adjustment.  A  test  network  designed  and  executed  for  another  reason  is  used 
here  to  demonstrate  the  relative  accuracies  attained  using  single  and  double 
runs  without  and  with  least-squares  adjustments. 

The  first  part  of  this  paper  shall  be  concerned  with  the  least-squares 
error  model,  while  the  remainder  shall  discuss  the  data  set  and  survey 
results. 

2.  LEAST-SQUARES  PROCEDURE 


The  adjustment  of  elevations  will  be  discussed  throughout  this  procedure, 
but  is  applicable  to  the  reduction  of  position,  deflection,  and  gravity 
components  as  well.  The  so-called  "SMOOTHED  ELEVATIONS"  are  meaned  to  obtain 
an  approximate  (preliminary)  elevation  for  each  station.  These  values  are  to 
be  used  in  the  computation  of  L  in  the  observation  equation: 

<*1  -  xj  - L  *  V*  (1> 

p  is  the  observation  weight 

is  the  unknown  value  of  the  i-th  station 

Xj  is  the  unknown  value  of  the  J-th  station 

L  is  the  difference  between  the  mean  ( preliminary)  values  and  the 
observed  values* 

or  L  -  [(J^  -  -  Oj))  (2) 

Constraint  Observation  Equation: 

(XA  -  L  «  V1)p 
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and  denoting  weights  'by  p^,  the  following  matrices  and 
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In  this  notation  the  observation  equations  are  given  by  AX=L. 

The  solution  in  the  sense  of  least  squares  is  X^A'pA)*1  (A'pL)*  (4) 

Given  two  legs  of  observations  from  the  IPS 
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Observations  equations  are  formed  thusly: 

To  express  it  in  equation  notation  and  eliminating  the  zero  elements 
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To  express  it  in  matrix  form 
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Results  of  least-squares  adjustment  showing  the  adjusted  elevations,  error 
terms  and  residuals* 
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3.  SAMPLE  DATA  SET 

In  1980,  the  U.S.  Army  Engineer  Topographic  Laboratories  observed  a  test 
network  at  the  White  Sands  Missile  Range  (WSMR)  as  a  cooperative  effort  with 
the  Defense  Mapping  Agency.  The  network  and  techniques  are  described  in 
detail  in  {Todd,  1981a  and  1981b),  The  unique  feature  of  this  data  set  is 
that  it  was  observed  between  known  survey  stations  in  something  approaching  a 
network  fashion,  which  permits  a  reduction  by  least-squares  methods,  and  gives 
a  reference  for  comparison  against  very  well  established  values. 

The  network  itself  is  shown  in"  Figure  1.  There  are  11  separate  test 
lines  containing  a  total  of  79  unique  survey  points.  Latitude,  longitude  and 
elevation  are  known  at  all  79  stations,  gravity  at  78  and  astrogeodetic 
deflection  of  the  vertical  at  47  of  the  stations.  Eighteen  of  the  79  stations 
were  treated  as  known  points  and  the  remaining  6l  stations  as  new  survey 
points.  Thirty  of  the  stations  are  at  the  Junction  of  two  or  more  traverses. 
Each  of  the  11  test  lines  was  run  once  each  way  (direct  and  reverse)  with  the 
system  in  a  vehicle  and  once  each  way  in  a  helicopter.  Only  the  vehicle  data 
set  is  used  in  this  illustration  of  least-squares  results. 

The  typical  time  between  ZUPTs  (zero- velocity  updates)  was  3. 5  minutes, 
and  the  average  length  of  a  one-way  traverse  between  updates  was  2.l8  hours. 

The  ETL  inertial  system  used  to  perform  the  observations  is  identical  to 
the  Defense  Mapping  Agency  IPS-1.  The  software  used,  however,  is  an  IPS 
program  modified  specifically  to  improve  the  accuracy  of  the  interpolation  of 
the  deflection  of  the  vertical.  That  on-board  software  is  referred  to  as  RGSS 
for  Rapid  Geodetic  Survey  System.  Position  and  gravity  accuracy  of  the  RGSS 
software  is  the  same  as  that  for  the  standard  IPS  program,  while  RGSS 
elevations  are  expected  to  be  slightly  inferior  unless  post-mission  processing 
is  employed. 

4.  PRELIMINARY  COMPUTATIONS 

The  RGSS  program  adjusts  the  latitude,  longitude  and  elevation  of  each 
intermediate  survey  station  between  the  end,  update  points  of  each  individual 
traverse.  The  adjustment  is  a  distribution  of  the  error  of  closure  of  each  of 


these  survey  parameters  at  the  terminal  update  station.  In  that  on-line 
smoothing  the  known  values  of  the  end  points  are  held  fixed. 

The  on-line  program  yields  output  from  which  the  anomalous  gravity  vector 
can  be  computed.  The  gravity  anomaly  is  given  as  the  Kalman  Z  accelerometer 
bias  DZ.  The  deflection  of  the  vertical  components  are  normally  computed  for 
the  IPS  software  by: 

q  ■  -(SUME  +  DN/4.85)  (6) 

=  -(SUMN  -  DE/4.85)  (7) 

where,  SUMN  and  SUME  are  the  Kalman  torque  about  north  and  east, 
respectively,  and  DN  and  DE  are  Kalman  north  and  east 
accelerometer  biases. 

However,  in  using  the  RGSS  software  the  deflections  are  computed  by: 

=  -DN/4.85  (8) 

=  DE/4.85  (9) 

The  relative  gravity  anomaly  and  deflection  of  the  vertical  must  be 
corrected  to  a  datum  based  on  a  known  value.  In  this  case  the  known  station 
was  the  beginning  point  of  the  traverse.  The  datum  shift  is  given  by: 

Fi  *  t±  +  C  (10) 

Where,  F^  =  the  observed  value  corrected  to  the  datum. 

f^  *  the  observed  value  at  the  ith  station. 

C  =  known  value ^  -  f 

The  gravity  anomalies  and  deflection  of  the  vertical  components  are  then 
adjusted  as  a  function  of  time: 


where  a  the  preliminary  value  of  the  gravity  anomaly  or  deflection 
component  for  the  least-squares  adjustment, 

Dj  *  the  error  of  closure  on  the  terminal  station  for  the  J  th  traverse. 

t^  »  the  time  of  observation  of  the  initial  update  station 

t ^  *  the  time  of  observation  of  the  ith  station 

t  ■  the  time  of  observation  of  the  terminal  update  station 

F^  *  the  observation  at  the  ith  station,  corrected  to  the  datum  in 
question  based  on  the  known  value  at  the  initial  update  station. 


The  statistical  term  used  in  all  computations  is  the  standard  deviation, 
one  sigma  level,  as  given  fcy  the  equation: 


where,  v  =  the  residual  formed  tor  the  difference  between  the  observed  and 

known  values. 

n  *  the  number  of  observations 
k  *  the  number  of  unknowns 

5.  THE  ADJUSTMENTS 

The  general  least  squares  adjustment  derxrlbed  in  this  paper  has  been  in 
use  for  some  time.  There  are  two  basic  options  available;  a  position 
adjustment  or  an  adjustment  of  elevation*  It  is  routine  procedure  to  store 
the  position  and  elevation  data  from  an  IPS  survey  in  a  permanent  automated 
data  file.  The  extraction  and  use  of  the  anomalous  gravity  vector  data  was  a 
special  procedure  for  this  test  project.  The  elevation  adjustment  option  was 
used  for  separate  adjustments  of  the  gravity  anomalies  and  the  deflection  of 
the  vertical,  components  with  only  slight  modifications  to  handle  the  formats 
for  input  and  output. 

The  data  sets  were  established  as  follows: 

Single  Run  -  All  residuals  from  22  runs  were  pooled.  The  residuals  were 
formed  by  differencing  the  individual  observed  and  smoothed  values  and  the 
known  values. 

Single  Run  Least-Squares  -  The  same  observations  used  above  were 
introduced  into  a  least-squares  adjustment  as  preliminary  values.  The 
residuals  were  formed  by  differencing  the  adjusted  and  known  values. 

Double  Run  -  The  observations  from  the  22  runs  were  sorted  into  11  double 
run  traverses.  The  smoothed  direct  and  reverse  observations  were  meaned  to 
yield  a  single  value  at  each  station  for  each  traverse.  The  residuals  were 
formed  by  differencing  these  smoothed  mean  values  and  the  known  values. 

Double  Run  Least-Squares  -  The  double  run  mean  values  were  introduced 
into  the  least-squares  adjustment  as  preliminary  values.  The  adjusted  and 
known  values  were  differenced  for  the  residuals. 

In  the  least-squares  adjustment  the  l8  known  stations  were  estimated  to 
have  a  standard  deviation  of  ±0.08  meters  in  latitude  and  longitude,  ±0.03 
meters  in  elevation,  ±0.3  seconds  of  arc  for  each  deflection  component,  and 
±0.1  mgal  for  the  gravity  anomaly. 

6.  RESULTS 

The  overall  results  from  the  analysis  of  the  U  data  sets  previously 
described  are  given  in  Table  1.  In  order  of  improved  accuracy,  the  data  sets 
are  single  run,  double  run,  single  run  with  least  squares  adjustment,  and 
double  run  with  least  squares  adjustment. 


(13) 


The  efficiency  quotient,  EQ,  is  computed  for  each  method  by: 

EQ  =  I/L 

where,  I  =  improvement  compared  to  single  runs 

L  =  level  of  effort  compared  to  single  runs 

The  level  of  effort  required  for  each  technique  varies  considerably 
depending  on  many  factors  such  as  size  of  project,  availability  of  good 
control,  spacing  of  stations,  etc.  On  a  project  of  100  or  more  stations,  the 
levels  of  effort,  compared  to  single  runs,  are  about  2  to  1,  1.2  to  1  and  2.4 
to  1  for  double  run,  single  run  network,  and  double  run  network  respectively, 
where  20%  of  the  stations  are  cross-connected  in  the  network. 

7.  SUMMARY 


While  post-mission  adjustments  of  inertial  positioning  data  have  been 
done  for  some  time,  several  problems  have  been  recognized.  Most  people  agree 
that  we  should  be  using  "raw"  inertial  data,  or  at  least  "raw"  survey  data  in 
the  adjustments,  but  rigorous  error  modeling  and  adjustment  techniques  for 
these  data  have  not  been  available.  So  we  now  adjust  "smoothed"  survey  data 
using  general  least-squares  adjustment  programs. 

There  has  been  a  general  acceptance  of  the  fact  that  the  inertial  survey 
systems  can  produce  sub-meter  position  and  elevation  accuracies  relative  to 
the  local  control,  but  the  anomalous  gravity  vector  has  been  generally 
neglected  except  by  the  U.S.  Army  Engineer  Topographic  Laboratories  with 
encouragement  from  the  Defense  Mapping  Agency. 

There  has  not  been  a  test  data  network  available  for  analysis  where  good 
survey  control  was  available  at  all  stations.  It  has  been  theorized  that 
post-mission  adjustment  must  improve  the  quality  of  the  network  surveys,  but 
the  magnitude  of  that  improvement  has  not  been  measured. 

The  Joint  DMA/ETL  effort  to  collect  and  analyze  data  from  a  pseudo 
network  observed  in  a  high  quality  test  area  will  not  solve  the  problem  of 
lack  of  rigorous  error  modeling  and  adjustment  techniques.  Efforts  to  produce 
such  rigorous  techniques  specifically  for  inertial  systems  must  be  given  more 
emphasis  and  encouragement. 

The  test  data  do  show  the  quality  of  survey  product  which  can  be 
obtained,  including  amazingly  good  recovery  of  the  anomalous  gravity  vector. 
This  aspect  has  been  ably  reported  by  ETL. 

The  purpose  of  this  paper  is  to  prove  that  the  post-mission  adjustments 
do  very  significantly  improve  the  accuracy  of  the  inertial  network  data.  This 
proof  is  given  in  the  table  of  results.  It  is  not  claimed  that  these 
techniques  can  remedy  the  problems  of  poor  quality  basic  survey  control, 
although  if  proper  relative  weighting  can  be  determined  then  some  improvement 
would  be  expected.  With  good  basic  control,  one  can  expect  about  a  25% 
improvement  in  overall  accuracy  of  double  run  data,  with  twice  the  survey 
effort,  compared  to  single  run.  That  yields  an  effectiveness  quotient  of 
0.12.5*  On  the  other  hand,  a  34%  improvement  in  overall  accuracy  is  to  be 
achieved  by  observing  a  single  run  network,  at  a  level  of  effort  1.2  times  the 
single  run  simple  traverses,  followed  by  post-mission  least-squares 
adjustment.  This  effectiveness  quotient  of  0.28  is  clearly  superior  to  the 


simple  double  run  method  in  those  cases  where  both  techniques  are  adaptable  to 
the  survey  situation. 

The  best  accuracy  is  obtained  by  the  double  run,  least-squares  adjusted 
network  data,  with  an  overall  kh%  improvement  in  accuracy.  The  level  of 
effort  is  2.U  times  that  required  for  simple  traverses  of  single  runs.  The 
effectiveness  quotient  of  0. 18  is  better  than  that  of  the  simple  double  run 
data,  but  not  as  good  as  the  0.28  of  the  single  run,  least-squares  adjusted 
data. 

For  those  relatively  large  projects  which  are  adaptable  to  network  type 
observations,  the  most  cost  effective  survey  method  is  single  run,  inter¬ 
connected  traverses  with  post-mission  adjustments. 
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